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ABSTRACT 

The  surface  ship  rolling  motion  equation  is  modeled  as  a  second  order  system,  with 
a  natural  frequency  of  con  —  0.4/.sec  and  a  dimcnsionless  damping  ratio  of  £  =  0.08.  The 
model  is  subjected  to  a  random  forcing  function,  which  has  a  Gaussian  probability  dis- 
tribution and  can  be  considered  as  "white  noise",  and  placed  into  State-Space  form.  State 
variable  feedback  of  roll  rate  is  applied  and  the  system  discretized  to  match  digital  con- 
trol. Roll  angle  time  histories  are  developed  for  a  range  of  feedback  gains  and  compared. 
Additionally,  steering  machinery  dynamics  are  modeled  by  a  first  order  system  and  time 
constants  varied  to  determine  the  effects  of  rudder  dynamics  on  the  feedback  system. 
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NOMENCLATURE 

A  Plant  Matrix 

bir        Element  B(2)  in  control  matrix,  sec-1 

B  Control  matrix,  [0  bir~]r 

k  Feedback  gain  constant,  seconds 

RMS    Root  Mean  Squared,  £-5(02) 

RRS   Rudder  Roll  Stabilization 

t         Time,  seconds 

tn        Time  at  "n-th"  time  sample,  (n)  Ts 

Ts       Sampling  time,  seconds 

u(t)    External  input,  continuous  form 

u(n)    External  input,  discrete  form 

X        State  variables,  [</>  0]r 

Sr        Rudder  control  angle 

C         Dimensionless  damping  ratio 

(f>        Roll  angle,  radians 

t         Steering  machinery  time  constant,  seconds 

(o„       LTndamped  natural  frequency,  radians/second 
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I.     INTRODUCTION 

Digitally  controlled  Rudder  Roll  Stabilization(RRS)  systems  are  currently  installed 
on  the  Hamilton  class  Coast  Guard  cutters  and  a  prototype  version  has  also  been  tested 
aboard  a  U.S.  Navy  Spruance  class  destroyer.  Observations  of  these  systems  have  shown 
that  up  to  a  40  percent  reduction  in  the  vessel's  root  mean  squared(RMS)  roll  amplitude 
value  is  obtainable.  Unfortunately,  the  maximum  benefits  achievable  by  RRS  are  not 
known.  This  deficiency  stems  from  the  performance  of  the  ship's  rudder  control  system, 
which  has  a  large  degree  of  non-linearity  and  also  operates  with  a  dead  band  at  small 
input  command  magnitudes  [Ref  1:  p.  46].  It  is  also  known  that  the  Spruance  and 
Ticondcroga  class  of  ships  have  a  steering  system  which  behaves  with  an  effective  time 
constant  that  varies  almost  inversely  with  amplitude  demand.  At  low  amplitudes,  2-3 
degrees,  the  system  can  be  modeled  as  having  a  3  second  time  constant  while  at  larger 
amplitudes,  7-9  degrees,  the  time  constant  reflecting  the  average  behavior  is  closer  to  1 
second.  In  addition,  the  four  separate  steering  systems  on  the  ship,  when  performing  at 
manufacturer's  specifications,  will  still  have  distinctly  different  dynamics  that  will  vary 
with  time  as  machinery  wear  and  crew  servicing  occurs. 

Additional  research  in  this  area  has  been  performed  by  Van  Der  Klugt  [Ref.  2]  and 
Powell  [Ref.  3],  to  name  a  few.  Van  Der  Klugt  focused  his  attention  on  the  incorpo- 
ration of  an  RRS  system  into  an  all  weather  autopilot  for  a  surface  vessel.  He  revealed 
that  the  RRS  system  needed  to  be  one  that  could  be  easily  modified,  to  meet  the 
changing  weather  conditions,  for  optimal  performance.  He  also  developed  guidelines  for 
the  proper  rudder  slew  rate  that  should  be  used  in  an  RRS  system  to  provide  maximum 
performance. 

Powell's  work  was  an  appraisal  of  RRS  and  compared  this  system  to  the  fin 
stabilizer  system.  What  this  report  exposed  was  the  need  for  a  larger  rudder  slew  rate 
and  increased  rudder  span.  It  also  discredited  proposed  ship  design  changes,  to  improve 
RRS  effectiveness,  such  as;  altering  the  rudder  Aspect  Ratio,  increasing  the  rudder  foil 
area,  angling  the  rudders  outboard  and  increasing  the  rudder  outreach.  He  also  discussed 
the  importance  of  developing  an  improved  steering  machinery  system  to  handle  the  RRS 
demands  and  agreed,  in  principle,  with  Van  Der  Klugt  on  the  rudder  slew  rate  limit. 

The  purpose  of  this  research  was  to  estimate  the  influence  of  steering  system  dy- 
namics, particularly  system  time  constants,  upon  the  RMS  roll  amplitude  of  a  ship  when 


the  rudder  is  used  in  a  stability  augmentation  system;  i.e.,  RRS  using  state  variable 
feedback  of  the  roll  rate.  While  most  research  in  this  field  concentrates  on  ship  design 
changes  and  rudder  slew  rate,  this  thesis  examines  machinery  response  to  the  control 
input  and  it's  effect  on  the  system  performance.  A  secondary  function  was  to  also  expose 
the  importance  of  an  adaptive  filter  which,  when  used  in  conjunction  with  the  digital 
RRS,  would  bring  uniformity  to  a  wide  variation  of  steering  system  dynamics  thus  pro- 
viding a  full  realization  of  the  stability  augmentation  capability  of  an  RRS  system. 
The  intended  procedure  to  accomplish  these  tasks  is: 

1.  Develop  a  system  math  model. 

2.  Discretize  the  model  to  match  digital  control. 

3.  Subject  the  model  to  a  generated  random  forcing  function  and  develop  a  roll  angle 
time  history. 

4.  Use  state  variable  feedback  of  roll  rate,  RRS,  and  compare  roll  angle  time  histories. 

5.  Develop  a  steering  system  model  and  incorporate  it  into  the  system  model. 

6.  Use  the  random  function  on  the  new  system  and  determine,  by  roll  angle  history 
comparison,  the  effects  q[  the  various  steering  system  dynamics  on  the  RRS  sys- 
tem. 


II.     VESSEL  DYNAMIC  MODELING 

As  a  surface  ship  travels  through  the  sea  it  is  subjected  to  three  primary  disturb- 
ances; the  wind,  waves  and  its  own  rudder  movement.  These  inputs  create  hydrodynamic 
forces  and  moments  that  cause  the  vessel  to  move  about  its  horizontal  axis,  or  roll.  The 
equation  that  best  describes  this  motion  is  a  linearized  fourth  order  state  equation  in- 
volving; sway  velocity,  yaw  rate,  roll  rate  and  roll  angle,  subjected  to  the  inputs  of  rudder 
deflection  and  sea  state.  Each  of  the  aforementioned  components  in  this  equation  has 
with  it  a  hydrodynamic  coefficient  that  is  variable.  The  determination  of  these  coeffi- 
cients can  be  quite  difficult  and  in  this  instance  the  sea  state  influence  coefficients  could 
never  be  found. 

Even  though  the  appropriate  coefficients  for  the  fourth  order  equation  could  not  be 
found,  it  would  still  be  possible  to  model  the  roll  dynamics  as  a  second  order  system. 
Granted,  this  model  would  be  less  stringent  then  the  fourth  order  equation  but  it  could 
still  provide  satisfactory'  and  realistic  results. 

A.     MODEL  OF  VESSEL  DYNAMICS 

A  surface  ship's  roll  dynamics  due  to  rudder  deflection  and  disturbances  can  be  ac- 
curately modeled  as  a  second  order  equation  of  motion  under  the  conditions  that  the 
model  provide  the  same  or  very  similar  results  as  the  actual  system  under  equal  condi- 
tions. Also,  the  model  results  cannot  be  considered  valid  when  used  during  instances 
where  the  actual  system  is  known  to  provide  conflicting  results.  As  an  example  consider 
that  a  second  order  model  may  be  stable  under  all  feedback  circumstances  while  the 
system  it  is  modeling  may  be  unstable  under  certain  conditions  of  state  variable  feed- 
back. Therefore,  the  model  results  cannot  be  used  under  the  conditions  of  feedback 
where  the  actual  system  is  unstable. 

The  second  order  model  for  the  equation  of  motion  describing  a  surface  ship's  be- 
havior in  roll  is: 

0  +  2Cco>  +  co2n<f)  =  b6r5r  +  u{t)  ( 1 ) 


where 


11 
(p  =  roll  acceleration,  d  4>\di 


4>  =  roll  rate,  d<\)\dt 

<f)  =  roll  angle 

£  =  dimcnsionless  damping  ratio 

con  =  natural  frequency 

b6r  =  rudder  scaling  factor 

Sr  =  rudder  deflection 

u(t)  =  input  disturbance 

The  variables  (  and  &>„  are  ship's  speed  dependent  and  were  selected  for  corre- 
spondence to  the  Spruance  class  destroyer  sea  trial  data  at  15  knots.  The  estimated 
magnitudes  were  £  =  0.0S  and  co„  =  0.4  /sec.  To  determine  the  magnitude  of  bir  one  must 
first  find  the  correlation  between  the  ship  roll  angle  and  rudder  deflection,  or  </>/<5f.  The 
sea  trial  data  was  again  used  and  a  value  of  0.15  roll  angle  (radians)  per  rudder  angle 
(radians)  was  found  from  steady  turn  observations. [Ref.  1:  pp.  45-46] 

Knowing  the  static  sensitivity  value,  0/(5,  ,  it  is  then  possible  to  return  to  the  second 
order  model  and  determine  bir  as  follows: 

<f>  +  2£g>„<£  +  con<f)  =  bdrSr 
at  steady  state  all  transients  decay  to  zero  and 

«^  =  b6r8r 
with 

—    * 

therefore  bir  =  0.024/  sec2. 

With  all  variable  coefficients  determined  and  the  disturbance  u(t)  acting  as  wind  and 
wave  action  the  second  order  model  becomes: 

if)  +  0.064c/>  +  0.\6<f)  =  0.024(5,  +  u(t) 


B.     STATE-SPACE  EQUATION 

The  second  order  model  can  be  converted  to  the  standard  State-Space  format: 


x  =  Ax  +  Bu 


by  recognizing  that: 


x  =  state  variable X-4>    ^T 

A  =  plant  matrix 

b  =  control  matrix 

u  =  noise  input 


The  State-Space  Equation  then  becomes: 


1 

">" 

+ 

"  0  " 

<5r  + 

"0" 

2C<V 

A 

_V 

f 

_1_ 

M(0 


(2) 


The  State-Space  form  provides  a  convenient  means  to  determine  and  analyze  the 
system's  response  to  various  manipulations.  The  addition  of  state-variable  feedback  or 
any  other  change  can  be  accomplished  with  relative  ease.  The  State-Space  representation 
can  be  quickly  converted  to  a  discrctized  form  to  accommodate  digital  control. [Ref.  4: 
p.  10] 

C.     DISCRETIZATION 

The  RRS  system  used  aboard  these  vessels  is  digitally  controlled  and  thus  our  model 
should  be  represented  and  used  in  this  manner.  The  governing  criteria  when  using  digital 
control  and  discrete  time  systems  is  the  sampling  frequency.  The  sampling  frequency  is 
how  often  the  control  system  senses  data  input,  or  the  amount  of  time  that  elapses  be- 
tween successive  data  samples.  The  choice  of  sampling  frequency  is  determined  by  the 
period  of  the  system  to  be  measured.  The  Nyquist  criteria  says  that  the  sampling  fre- 
quency must  be  no  less  than  two  times  the  maximum  frequency  of  the  system  being 
measured  [Ref.  5].  It  has  been  determined  that  the  period  of  oscillation  for  the  Spruance 
destroyer  at  15  knots  is  15.3  seconds  which  gives  it  a  frequency  of  0.0653/  sec  [Ref.  1: 
p.  44].  A  sampling  time  step  of  0.333  seconds  was  chosen,  which  is  greater  than  45  times 


the  frequency  of  the  system  and  allows  for  a  complete  reconstruction  of  the  continuous 
input  disturbance  function. 

When  dealing  with  discretization,  state  variables  are  predicted  for  each  successive 
time  step  by  the  previous  step's  state  variables  and  the  previous  input.  A  numerical 
representation  of  this  is: 


x(n  +  1)  =  <&(Ts)x(n)  +  r(Ts)u(n) 


where 


<&(Ts)  =  eATs 

r{Ts)  =  [eArs-\]A~xB 

I  =  identity  matrix 

A  =  plant  matrix 

Ts  =  time  step 

x  =  state  variables 

D  =  control  matrix 

u  =  input  disturbance 

n  =  integer  index,  n  =  0,1,2,3,.... 

By  using  this  type  of  approximation  a  recursive  relationship  is  developed  and  the 
state  variable's  time  history'  can  be  easily  calculated  by  a  computer  program. [Ref.  6] 

D.     FEEDBACK 

State  variable  feedback  can  be  used  to  alter  the  dynamic  response  of  a  system.  In 
this  case  the  desired  effect  is  to  decrease  the  roll  angle  magnitude  of  the  ship  by  feeding 
back  the  roll  rate  increased  by  a  multiplicative  gain. [Ref.  4:  p.  222] 

To  maintain  the  accuracy  of  the  second  order  model  it  is  necessary  to  look  more 
closely  at  the  actual  fourth  order  Equation  of  Motion  for  this  system  to  determine  if  and 
where  that  system  may  achieve  maximum  stability  or  become  unstable.  An  estimated 
root  locus  plot  of  the  fourth  order  system  has  been  developed  and  the  optimum  value 
or  position  of  greatest  stability  occurs  when  the  dimensionless  damping  ratio,  (,.  is  equal 


to  0.4  for  a  ship  speed  of  15  knots  [Rcf.  1:  p.  47].  This  implies  that  the  maximum 
damping  ratio  for  the  feedback  compensated  system  model  can  not  be  expected  to  ex- 
ceed C  —  0.4.  Figure  1  displays  the  estimated  root  locus  plot  for  the  ship. 
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Figure   1.      Root  Locus  plot 


The  State-Space  Equation  then  becomes: 

w 


-con       -2Ccon  ~kb6r\  10 


+ 


«(0 


(3) 


where 


-2C<rW„  =  -2Cco„  -  kb 


6r 


and  the  maximum  value  of  £e  =  0.4.  Therefore  the  maximum  value  of  the  gain  k  is  10.66. 
This  system  was  discrctized  as  discussed  previously. 

E.     RUDDER  CONTROL  MODEL 

The  ship's  electro-hydraulic  rudder  control  system  behaves  non-lincarly  and  is  also 
accompanied  by  a  dead  band  at  small  rudder  input  commands.  Since  the  dead  band  is 


very  small  relative  to  the  full  scale  it  will  be  ignored  in  the  model  and  the  non-linear  re- 
sponse will  be  characterized  by  a  first  order  system  with  the  various  time  constants  re- 
presenting the  average  first  order  behavior  for  different  amplitudes  of  the  machinery. 
The  first  order  equation  for  the  steering  system  is: 

5r  =  4-  (*«  -  *r)  (4) 

(5r  =  rudder  deflection  rate,  sec~ 

x  =  lime  constant,  seconds 

5rc  =  helm  command 

Sr  =  rudder  deflection 


From  shipboard  machinery'  performance  tests  it  was  determined  that  the  time  con- 
stants associated  with  the  steering  system  ranged  from  one  to  three  seconds  and  these 
values  will  be  used  in  this  discussion  [Rcf  1:  p.  46].  Considering  that  a  ship  usually 
travels  in  a  straight  path,  it  was  decided  to  use  a  zero  degree  rudder  command  from  the 
ship's  helm.  This  is  not  a  poor  assumption,  most  vessels  remain  on  the  same  course  for 
long  lengths  of  time.  This  reduces  the  rudder  model  to  the  time  constant  and  the  rudder 
deflection  and  deflection  rate.  The  actual  ship  steering  machinery  has  an  approximate 
rudder  demand  rate  limit  of  6  degrees  per  second  [Rcf.  1:  p.  47].  In  the  studies  described 
herein,  it  was  assumed  that  required  rudder  rates  did  not  involve  the  influence  of 
control-rate  limiting.  Figure  2  displays  the  rudder  deflection  with  various  time  constants 
to  a  unit  step  input. 
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Figure  2.      Rudder  Response 


F.     THE  COMPLETE  SYSTEM 

When  the  rudder  deflection  is  included  in  this  system,  to  account  for  the  rudder  dy- 
namics, it  assumes  the  role  of  a  third  state  variable  and  the  <&(Ts)  and  r(7s)  discussed 
previously  in  the  second  order  system  are  now  represented  in  a  third  order  equation.  The 
complete  system  including  feedback  and  the  steering  system  characteristics  when  placed 
in  State-Space  form  becomes: 


*1 

0 

1 

0 

,<K 

0 

*}- 

2 

-2Co>„ 

b6r 

v\+ 

1 

V 

0 

-A/t 

-1/t 

kJ 

0 

"(') 


(5) 


III.     RESULTS 

The  random  function  used  as  the  disturbance  from  wind  and  wave  action  was  de- 
veloped by  the  Monte  Carlo  method  [Rcf.  7].  The  characteristics  of  this  function  are 
that  it  has  a  Gaussian  probability  distribution  and,  because  it  represents  a  wide-band 
random  input  function,  can  be  considered  as  white  noise.  In  essence,  the  wide-band  in- 
put function  passes  through  a  band-pass  filter  (the  ship  math  model)  and  produces  a 
narrow-band  output  similar  to  that  of  the  actual  ship's  behavior.  One  must  recognize, 
however,  that  the  sea  is  not  considered  "white"  and  that  this  white  noise  function  was 
used  to  simplify  the  problem.  The  same  disturbance  was  used  on  all  the  systems  thereby 
validifying  the  comparisons.  A  graphical  representation  of  the  disturbance  function  is 
shown  in  Figure  3. 
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Figure  3.      Disturbance  function 
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A.     FEEDBACK  EFFECTS 

To  determine  the  efTects  of  roll  rate  feedback,  upon  the  system  with  ideal  steering 
characteristics,  a  scries  of  output  responses  were  developed  for  a  gain  magnitude  range 
of  zero  to  10.66.  This  corresponds  to  a  damping  ratio  spectrum  of  0.08(no  gain)  to  0.4(k 
=  10.66).  Figure  4  displays  the  roll  angle  RMS  for  this  spectrum  and  as  expected  the 
RMS  magnitude  significantly  decreases  from  3.641  to  1.692,  a  53  percent  reduction.  Also 
plotted  on  this  graph  is  the  expected  RMS  behavior  of  a  second  order  system  subjected 
to  white  noise,  RMS  oc  l/jc,  .  A  comparison  between  the  two  curves  shows  that  the 
model  results  are  slightly  higher.  This  can  be  accounted  for  by  round  off  error  in  both 
the  discretization  and  recursion  relationship  calculations.  Since  the  output  response  of 
the  ship  model  compares  very  closely  to  the  theoretical  response  of  the  second  order 
system  to  white  noise,  it  can  be  assumed  that  the  ship  model  is  valid. [Ref.  8:  pp.  503-505] 
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Figure  4.      RMS  vs.  Damping  ratio 


Figures  5  and  6  are  the  actual  roll  angle  time  histories  for  the  uncompensated  and 
fully  compensated  systems.  From  these  figures  it  is  obvious  that  not  only  has  the  RMS 
value  decreased  but  that  the  peak  magnitudes  have  been  drastically  reduced.  It  is  the 
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reduction  in  peak  magnitude  that  is  really  the  most  beneficial  aspect  of  RRS,  there  is 
always  going  to  be  some  rolling  motion  and  if  the  size  of  the  roll  angle  can  be  minimized 
then  the  ship  and  her  crew  can  carry  out  it's  mission  more  effectively. 
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Figure  5.      Uncompensated  Roll  Angle  History,  Spruance  class  Destroyer 
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B.     STEERING  SYSTEM  EFFECTS 

Time  constants  of  zero  to  three  seconds  were  chosen  for  use  with  the  system  in- 
volving the  maximum  gain  magnitude  of  10.66  (damping  ratio  of  0.4).  Figure  7  displays 
the  relationship  between  the  output  RMS  of  the  assumed  ship  model  to  the  individual 
time  constants.  The  RMS  magnitude  increases  with  the  increase  in  the  time  constant 
thereby  substantiating  an  intuitive  feeling  that  steering  machinery  dynamics  adversely 
effect  the  RRS  system  and  more  importantly,  that  the  influence  of  the  steering  machin- 
ery can  be  quantified. 
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By  taking  a  closer  look  at  equation  5  one  can  fully  realize  the  effects  of  the  ma- 
chinery time  constant  on  the  system  behavior.  Figure  8  is  a  plot  of  the  characteristic 
roots  of  that  equation  based  on  various  time  constants  and  from  this  plot  it  is  clear  that 
as  the  time  constant  increases  the  effective  system  damping  ratio  is  decreasing.  This  is 
more  clearly  displayed  in  Figure  9  where  the  damping  ratio  is  plotted  versus  the  time 
constants.  It  is  because  of  this  decrease  in  system  damping,  brought  on  by  the  addition 
of  the  time  constants,  that  increases  the  RMS  value  of  the  roll  angle. 
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Figure  8.      Characteristic  Roots  of  Equation  5 
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Figure  9.      Damping  Ratio  VS.  Time  Constants 

Figures  10  and  11  display  the  output  response  of  the  zero  and  1.5  second  time  con- 
stant systems  using  the  maximum  gain  of  k=  10.66.  From  these  figures  we  see  an  in- 
crease in  peak,  magnitude  and  roll  angle  RMS  as  the  time  constant  increases  with  no 
other  significant  differences.  Table  1  shows  the  percentage  of  RMS  reduction  for  time 
constants  ranging  from  1  to  3  seconds.  Without  steering  machinery  dynamics,  the  po- 
tential RMS  roll  reduction  due  to  the  assumed  state-  variable  feedback  was  53  percent. 
The  effect  of  a  one  second  time  constant  reduced  the  potential  roll  reduction  by  ap- 
proximately 6  percent  while  a  three  second  time  constant  reduced  the  potential  roll  re- 
duction by  24  percent.  The  conclusion  is  that  "sluggish"  steering  machinery  dynamics 
has  a  perceptible  effect  upon  ship  roll  stabilization  by  the  use  of  the  rudder. 
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Figure    10.       Roll  Angle  History,  Ideal  Steering  Machinery 
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Table  1.     RMS  AND  PERCENT  REDUCTION  VAL- 
UES FOR  TIME  CONSTANT'S 


Tunc  constant  in  seconds 

RMS 

Percent  re- 
duction 

0.0 

1.692 

53.5 

0.5 

1.719 

52.3 

1.0 

1.794 

50.7 

1.5 

1.897 

47.8 

2.0 

2.008 

44.8 

2.5 

2.113 

41.9 

3.0 

2.240 

38.4 

Figures  12  and  13  display  rudder  deflection  time  histories  for  the  fully  compensated 
system  with  two  different  time  constants.  The  point  to  be  made  is  that  as  the  time 
constant  increases  the  magnitude  of  the  rudder  deflection  decreases  because  the  system 
becomes  less  responsive.  Thus  the  rudder  becomes  less  effective  and  it  then  cannot  fully 
support  the  RRS  system. 
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Figure   12.       Rudder  Activity,  Machinery  Time  Constant  of  U.5  Seconds 
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IV.     CONCLUSIONS  AND  COMMENTS 

It  is  obvious  from  the  results  that  steering  machinery  dynamic  characteristics  ad- 
versely effect  the  Rudder  Roll  Stabilization  system.  Even  with  this  hindrance  the  system 
can  still  reduce  the  roll  angle  RMS  and  peak  magnitude  values  significantly. 

Up  to  this  point  no  mention  of  the  effects  of  the  RRS  system  upon  the  rudder  itself 
has  been  made.  Figures  12  and  13  show  that  the  rudder  is  moving  often  and  very  quickly. 
Concerns  directed  towards  machinery  wear  and  reliability  are  justified  and  focus  should 
be  placed  on  designing  these  systems  to  handle  the  increased  activity.  If  no  design  sol- 
ution can  be  found  then  guidelines  should  be  developed  suggesting  that  the  RRS  system 
only  be  operated  in  harsh  sea  states  thus  reducing  the  overall  operating  time  and  possi- 
bility of  material  failure. 

The  choice  of  the  maximum  system  gain  was  made  purely  on  the  basis  of  stability. 
Van  Dcr  Klugt  in  his  paper,  "Rudder  Roll  Stabilization:  The  Dutch  Solution",  states 
that  there  is  a  point,  irrespective  of  the  rudder  speed,  where  the  RRS  system  can  longer 
reduce  the  roll  activity.  This  point  is  a  rudder  speed  of  approximately  20  deg/sec  which 
corresponds  to  a  compensated  damping  ratio  of  0.52  in  this  system. [Ref.  2  :  p.  85]  In 
this  research  that  point  could  not  be  reached  because  the  gain  needed  to  get  there  would 
have  driven  the  fourth  order  equation  of  motion  unstable,  as  seen  by  Figure  1. 

On  the  other  hand,  Powell's  paper  ,  "Rudder  Roll  Stabilization:  A  Critical  Review," 
has  determined  that  a  rudder  speed  of  approximately  10  deg/sec,  which  is  an  average 
value  taken  from  the  results  of  various  sea  states  and  wave  encounter  angles,  is  where 
the  point  of  no  further  roll  reduction  takes  place  [Ref.  3:  pp.  2.250].  Using  that  criterion, 
one  sees  that  the  second  order  model  used  in  this  research  is  reasonable  and  produces 
qualitativly  correct  results. 

Another  conclusion  that  can  be  drawn  is  that  each  piece  of  machinery  is  going  to 
behave  differently.  This  fact  promotes  the  development  of  an  adaptive  filter  that  can  be 
used  to  bring  uniformity  to  the  wide  range  of  machinery  dynamics  thus  enabling  the 
RRS  to  reach  its  full  potential  as  a  stability  augmentation  device. 


APPENDIX      THE  INVERSE  PROBLEM 

The  original  idea  behind  this  research  was  to  use  the  roll  angle  time  histories  developed 
by  the  sea  trial  team  of  Baitis  and  Schmidt  ct.  al.  and,  working  in  an  inverse  manner, 
develop  the  sea  state  forcing  function  that  caused  this  motion.  This  idea  had  good  mer- 
its. The  roll  time  history  data  was  taken  when  the  ship  was  on  a  steady  course  with 
minimal  rudder  movement,  so  the  roll  data  could  reasonably  be  considered  from  a  rudder 
fixed  condition  which  would  mean  that  all  motion  was  induced  from  the  sea.  The  coef- 
ficients for  the  fourth  order  state  equation  were  known  with  good  confidence,  except 
those  for  the  sea  state,  so  it  seemed  that  the  inverse  problem  could  be  solved  by  working 
backwards  through  the  state  equation.  What  stood  in  the  way  was  the  fact  that  three 
out  of  the  four  state  variables  were  either  measured  inaccurately  or  not  measured  (ob- 
served). The  yaw  rate,  sway  velocity  and  roll  rate  were  missing  from  the  available  data 
base  and  without  them  it  seemed  that  any  invcrsing  process  would  be  doubtful. 

These  concerns  were  substantiated  by  Gao  and  Hess  who  describe  an  inverse  tech- 
nique which  is  an  iterative  process  requiring  that  all  output  states  were  measurable  and 
that  the  number  of  inputs  must  equal  or  exceed  the  number  of  outputs  [Rcf.  9J.  Since 
neither  of  these  conditions  were  met,  the  research  steered  away  from  this  path. 

In  retrospect,  a  solution  to  this  problem  may  have  been  to  use  the  actual  sea  trial 
roll  data  in  conjunction  with  a  second-order  linear  ship  model  to  develop  an  approxi- 
mation to  a  forcing  function.  Then  the  input  function  could  be  re-applied  to  the 
second-order  model  system  with  state  variable  feedback  and  steering  system  time  con- 
stants, so  that  roll  time  histories  could  be  developed.  These  time  histories  could  then  be 
compared  and  the  impact  of  the  steering  system  dynamics  brought  out.  This  approach 
is  left  for  future  consideration  in  follow  on  studies. 
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